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Abstract 



We calculate the dynamical response of two-layer Josephson arrays to a 
bias current applied into only one layer (primary layer), in a perpendicu- 
lar magnetic field. The pancake vortices in the two layers sometimes form 
into flux lines which move as a unit under the influence of the bias current. 
The lines break apart into independent pancakes in sufficiently high currents, 
sufficiently weak interlayer coupling, or possibly sufficiently high fields. We 
discuss the relevance of these results to recent experiments in high-T c super- 
conductors in the flux transformer geometry. 
PACS numbers: 74.50.+r, 74.60.Ge, 74.60.Jg, 74.70.Mq 
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I. INTRODUCTION 



The so-called "flux-transformer" geometry provides a novel way to study vortex 

order in Type-II superconductors. In this technique, the external current density J is injected 
into and parallel to the upper ah plane of a high-T c superconductor. The voltage drop is 
measured across both the upper plane (the "primary voltage") and the lower plane (the 
"secondary voltage"). For applied magnetic field B||c, the geometry is sensitive to the 
vortex order, especially to the interlayer coupling between the so-called vortex "pancakes" 
in the 2D layers. Typically, this coupling will cause the pancakes to form into flux lines. 
The J x B Magnus force will exert a torque on such a flux line. If the interlayer coupling is 
strong enough, the flux line will move as a unit under this force, producing equal primary 
and secondary voltages. But at sufficiently large J, for any large anisotropy, the influence 
of the torque will exceed that of the interlayer coupling, and the flux line will break. Even 
at small applied currents, the flux line can sometimes be broken by thermal noise. In both 
cases, the line breaking shows up in the difference between the primary and the secondary 
voltages, which increases quickly from zero as the pancakes decouple.) This cross-over from 
3D flux lines to 2D vortex pancakes has attracted great experimental and theoretical interest 

M3- 



In this paper, rather than treating the actual high-T c superconductor directly, we model 
the T = flux line dynamics in the flux-transformer geometry, using a two-layer Joseph- 
son junction array. While far from a realistic high-T c material, this model is numerically 
tractable [ 12jl3l , and in addition could be directly studied experimentally. We find that in 
this flux-transformer geometry, the vortex lines do exhibit a 3D-2D decoupling transition at 
a critical vortex velocity, just as in the high-T c materials. We map out this transition as a 
function of anisotropy and applied current, and find behavior which qualitatively resembles 
experiment. Our calculations show that this transition can be caused by Josephson coupling 
alone, without other types of interactions (e. g., electromagnetic interactions |],[|, Coulomb 
interactions [0]) between the pancakes on different layers. 
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II. MODEL 



We consider a two-layer square arrangements of Josephson junctions (cf. Fig. 1). All the 
junctions, both within and between the planes, are assumed to be overdamped, resistively 
shunted junctions (RSJ model). All intralayer and interlayer junctions are assumed identical, 
with critical currents and resistances I c , R and I C >,R respectively. The ratio a = I c >/I c 
represents the array anisotropy, and is the analog of the high-T c mass anisotropy. 

To simulate the flux-transformer geometry, we inject a bias current along the y direction 
into only the upper layer (the "primary" layer; cf. Fig. 1) ||. We denote the other the 
secondary layer. To introduce vortices, we apply a magnetic field B||z, i. e., perpendicular 
to the layers. The Magnus (J x B) force from the transport current will drive these vortices 
in the x direction. 

The equations of motion for this system take the standard form i^- = Vij/Rij + 
J c;ii sin(0j - 4>j - Aij, Vij = VI - Vj = (h/2e)(d/dt)(<f>i - 0,), and = h,ext- Here 

Ijj is the total current from grain % to grain j; Rij and I c; ij are the shunt resistance and 
critical current of the Josephson coupling between grain i and grain j; 0« is the phase of 
the order parameter on grain i; and Ii ]ex t is the external current fed into grain i. and 

are the voltage difference and magnetic gauge phase factor between grain i and grain j. 

= (27t/$ ) A • dl, where $ = hc/(2e) is the flux quantum, A is the vector poten- 
tial for the applied magnetic field and Xi is the position of the center of grain i. The first 
equation describes the total current between grains i and j as the sum of two terms: an 
Ohmic current through the shunt resistance, and a Josephson current. The second is Joseph- 
son equation relating voltage and phase, while the last is just Kirchhoff 's Law, describing 
current conservation at the i th node. 

As boundary conditions, we introduce a uniform current I^ext = I into each boundary 
grain in the primary layer and extract it from the other edge (cf. Fig. 1). In the secondary 
layer, we use free boundary conditions in the y direction. In both layers, periodic boundary 
conditions are adopted in the x direction. Am = for an interplanar junction, while for an 
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intraplanar junction Y^piaquette — ^"ff = 27r/ (S is the plaquette area) [|T^]. We consider 
only / = 1/N 2 so that there is only one vortex (i. e. one flux quantum) piercing each layer. 
We denote these the primary and secondary vortices. 

We solve the coupled equations of motion as described previously Jl5[], using a fourth- 
order Runge-Kutta algorithm. During the calculation of the dynamical state {<fii(t)}, we 
keep track of the vortex motion in the array, in terms of the vortex number n of each 
plaquette. n is defined by Yypiaquettei'Pi ~ 4>j ~ Aij) + — 2vrn, where (0j — 4>j — Aij) is in 
the range (— tt, it] and the summation is taken in the counterclockwise direction. 



III. I - (V) CHARACTERISTICS 

Our / — (V) characteristics typically fall into one of five regimes, (i) At the lowest 
currents, (V) = 0: both primary and secondary vortices are pinned by the periodic potential 
of the array Jl6].[17] • (ii) At slightly higher currents, in some cases, the primary vortices move, 
while the secondary vortices are still pinned, (iii) In other arrays, at such currents, both 
primary and secondary vortices are depinned and move together as a single flux line, (iv) 
At still higher currents, in some arrays, both vortices are depinned, but move independently, 
as separate "pancakes", (v) For / > J c , the entire Josephson lattice is depinned, and the 
/ — (V) characteristics are dominated by single-junction effects. 

Fig. 2 shows the calculated I — (V) characteristics of a 2 x 10 x 10 (N = 10) array at 
/ = 0.01 and three values of a. In the flux-flow regime [cases (ii)-(iv)], the J— (V) curves are 
dominated by vortex motion. In all cases, whether moving or fixed, these vortices remain 
in the center of the array (i. e. the 5 th row), presumably the region of lowest energy. (In our 
19 x 19 array calculations, they stay in the 10 th row.) 

The critical current for the onset of voltage is very sensitive to the anisotropy a. When 
a = 0, Fig. 2(b) shows that the primary vortex is depinned at the bias current I cp ~ 0.1/ c , 
while the secondary vortex remains fixed up to a bias current I cs w 0.28/ c . At a = 0.01, 
these depinning currents are respectively 0.16J C and 0.26/ c . But when a = 0.04, the primary 
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and the secondary vortices are depinned simultaneously near 0.17J C , and move thereafter as 
a single flux line in the array. This line, however, breaks apart near = 0.36/ c , above 
which the primary vortex moves faster than the secondary vortex. For all three a values, 
the individual junctions in the primary layer go normal near J CJ = 1.06/ c . 

In the zero-coupling limit, the value I cp = 0.1I C [cf. Fig. 2(b)] is the depinning current 
of a single vortex in a 2D array fl6| . This can be understood by the following argument. 
Below Icp, there is zero voltage drop across every junction. Hence, no current passes through 
the resistive shunts between layers. Moreover, since a = 0, there is no interlayer Josephson 
current. Therefore, all the bias current passes only through the primary layer, and I cp 
must be the vortex depinning current for a single-layer array. Once the primary vortex has 
been depinned, there is some voltage drop between the layers. This voltage drop produces 
currents into the secondary layer via the the interlayer shunt resistances. For sufficiently 
large /, these currents generate a Magnus force which depins the secondary vortex. 

At sufficiently small nonzero a, we still find that I cp < I cs . But above a critical a c , 
Icp — Ics, and the two vortices move together as a single flux line for I > I cp . a c seems also 
to be size-dependent. 

We can crudely understand all of these phenomena from a simple argument. The idea is 
illustrated schematically in Fig. 3, which shows a side view of an N = 11 array with transport 
current injected from the left. The primary and secondary vortices - both located in the 
middle row of the array - are indicated by crosses. Since the vertical junctions have critical 
current al c , the total current passing through these junctions into the secondary layer, 
under zero voltage conditions, cannot exceed ~ \aNI c to the left of the vortex. Assume 
that the transport current injected into the primary layer will divide equally between the 
two layers until this limit is exceeded. This occurs at a transport current I eq = aNI c . From 
this picture, we can estimate the conditions on iV and a such that the two vortices will 
be depinned simultaneously. Namely, if Na/2 > 0.1, then both vortices will be depinned 
simultaneously at current below this limiting current and move together thereafter under 
the bias current drive. Conversely, if Na/2 < 0.1, the vortices will be depinned at different 
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currents, and will move more or less independently. (Note that if there were no resistive shunt 
coupling between layers, then at a = value, the two layers would be entirely independent 
the secondary vortex would never be depinned.) 

For iV = 10, the above argument gives 0.02 < a c < 0.03, consistent with our numerical 
results. For an iV = 19 array, the corresponding value of a c , both analytically and numeri- 
cally, is in the range 0.01 < a c < 0.02. To test this picture further, we have also calculated 
the current distribution in the array for a variety of parameters. The distribution of currents 
in the vertical junctions indeed conforms with the estimates given above. 

This picture may possibly explain the breaking of the flux line at large enough /. As / 
increases, the line moves faster and faster, but the secondary vortex will lag the primary one 
by an increasing distance. Eventually, this distance may be such that (in a small array) the 
secondary vortex will start to feel the influence of the images of the primary vortex. Above 
this current, the primary and secondary vortices should break apart and move separately. 

These arguments suggest two conclusions. First, for large enough N, the primary and 
the secondary vortex will always form a single flux line - that is, will be depinned and move 
together - for an arbitrarily weak but nonzero interlayer Josephson coupling a. Secondly, 
in a sufficiently large array, the line should remain unbroken all the way up to the single- 



junction critical current I c ||18|| . We have carried out a few calculations for larger arrays, 
which do support these conclusions. 

We can also speculate how the 3D/2D transition would depend on magnetic field strength. 
Since / = 1/N 2 in our simulations, decreasing the array size may be comparable to decreas- 
ing the field strength. Our simulations at different N, combined with the above arguments, 
suggest that a c increases with increasing / - that is, that flux lines are more easily broken 
in high flux density superconductors. Hence, for a given weak interlayer coupling, there 
may be an upper limiting field, above which there are no 3D flux lines at all, but only 2D 
pancake vortices. Such an upper limit may have been observed in experiments J7| and has 
been predicted theoretically [0,Bf. 
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IV. TIME-DEPENDENT BEHAVIOR 



Fig. 4(a) shows both V p (t) and V s (t) for N = 10, I = 0.2I C , and a = 0. For this case, we 
find that the primary vortex moves across the array in the 5 th row, while the secondary vortex 
remains fixed in the plaquette (10, 5). Correspondingly, (V p ) ^ but (V s ) = 0. (Similar 
results are found for larger arrays.) The peaks in V p (t) correspond to the translation of 
the primary vortex by one plaquette. The corresponding peaks in V s (t), resulting from the 
weak resistive coupling between the layers, correspond to oscillations of the secondary vortex 
around an equilibrium position. Since the interlayer resistive coupling is very weak, there is 
a constant time interval between consecutive peaks in V p (t), implying that the velocity of the 
primary vortex does not depend on its proximity to the secondary vortex, and suggesting 
that the intervortex interaction nearly vanishes under these conditions. The time interval T 
between peaks satisfies the Josephson relation (V) — Ye~Wf-> wnere (Vp) is the time-averaged 
voltage difference across the primary layer. This is consistent with the picture of a vortex 
circling the array once in an interval NT. 

Fig. 4(b) shows the corresponding voltages at N — 10, / = 0.16/ c , and a = 0.01. In 
this case, the primary vortex is depinned, traveling in the 5 th row of the array, while the 
secondary vortex remains fixed in plaquette (10, 5), as in the a = case. Once again, both 
V p (t) and V s (t) exhibit many peaks, with the same interpretation as in 4(a). But now the 
interval between peaks on V p (t) is variable, implying a plaquette-dependent primary vortex 
velocity. The interval is smallest, and the average vortex velocity greatest, when primary and 
secondary vortices interval are directly superposed, and the mutual attraction is strongest. 
Because of the periodic boundary conditions, V p (t) is periodic, with N peaks in each period. 
The peak structure in V p (t) proves that an attractive interaction between vortices can be 
produced by interlayer Josephson coupling alone. 

Fig. 5 shows V(t) and the vortex paths for N = 10, / = 0.01, a = 0.04, and / = 0.18/ c . 
In this case, the two pancake vortices form a single flux line which moves as a unit at low 
bias current. Both V p (t) and V s (t) are characterized by periodic double-peaked structure 
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[cf. Fig. 5(a)], which can be understood from the vortex spatial path shown in Fig. 5(b). 
The larger peaks of V p (t) and the smaller peaks of V s (t) correspond to the translation of 
the primary vortex by one plaquette, while the other peaks arise from the translation of the 
secondary vortex. V p (t) and V s (t) have the same frequency. Although the two vortices move 
as a unit, the vortex line is at an angle to the c axis: the primary vortex leads the secondary 
vortex on average by about half a plaquette. In simulations at iV = 19, we have found 
that the primary vortex can lead the secondary vortex by as much as two plaquettes. As a 
increases, this horizontal distance decreases, showing that the interlayer vortex interaction 
is increasing. 

Note also that V p (t) and V s (t) have slightly different waveforms, although they have the 
same frequency. This means that the flux line formed by the two pancakes is not rigid. Our 
picture is that the flux line moves as a unit across the array in a periodic motion, but that 
the two constituent pancakes may oscillate around their steady state positions on the flux 
line. The period of such an oscillation would equal that of the flux line motion. 

At / = 0.36/ c , for this coupling, the flux line breaks into independent primary and 
secondary pancakes. Above this current, (V p ) > (V s ), indicating that the primary vortex 
moves faster than the secondary one. Such a critical vortex velocity for flux line breaking, 
i. e., a 3D to 2D crossover, has been reported in high-T c superconductors [fLl]| . 

V. DISCUSSION AND CONCLUSIONS 

We have studied the dynamical response of two-layer Josephson arrays to a bias current 
applied in the flux transformer geometry. Although this model is far from a realistic high-T c 
superconductor, our numerical results agree qualitatively with some experimental studies 
on such materials. Thus, they may be useful in constructing more realistic models for the 
dynamics of high-T c materials in this geometry. 

Our results show that interlayer Josephson coupling alone is enough to bind pancake 
vortices into vortex lines. These lines move as a unit under the drive of an external bias 
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current, but may break up into 2D pancake vortices and move independently at high enough 
currents, provided the magnetic field is not too small. Our calculations also suggest that, 
for any given current, there may exist an upper magnetic field limit above which 3D flux 
lines are unstable in the flux transformer geometry. 

Our results suggest several areas for future work. First, it would be of interest to extend 
this work to larger magnetic fields and finite temperatures (included via Langevin noise), in 
order to study the 3D-2D transition at finite fields and temperatures. Likewise, the inclusion 
of a finite inductive coupling between plaquettes [0,0 might lead to novel behavior not 
seen in the present, low-screening limit. Finally, the predicted effects might be sought 
experimentally, not only in high-T c materials, but directly in two-layer Josephson arrays 
which could be made in the geometries described here. 
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FIGURE CAPTIONS 



1. Schematic diagram of a 2 x 5 x 5 (N = 5) two-layer Josephson junction array. Each 
intersection denotes a superconducting grain. Each grain is coupled to its four in- 
tralayer neighbors, and to the nearest grain in the adjacent layer, by a resistively 
shunted Josephson junction (RSJ). An external magnetic field B is applied perpendic- 
ular to the layers in the z direction. A bias current I is injected into each grain at one 
edge of the upper layer (the "primary layer"), and extracted at the opposite edge in 
the y direction. V p and V s denote the voltage drops across the primary and secondary 
layers (averaged over x), as shown. There are periodic boundary conditions in the x 
direction. 

2. I — (V) characteristics for an iV = 10 at / = 0.01, and three different values of a. In 
each case, the two curves denote V p and V s . (a) a = 0; (b) a — 0.01; (c) a = 0.04. 

3. Schematic view of the array, viewed along the x direction, and showing the current 
direction in junctions parallel to y axis in an N = 11 array. Arrows denote the direction 
of current flow; crosses denote the vortex pancake positions. 

4. V p (t) (full curves) and V s (t) (dotted curves) in N = 10 arrays at / = 0.01 for two 
coupling ratios and bias currents as indicated. 

5. Time-dependent voltage traces and vortex path for an N = 10 two-layer array at 
/ = 0.01, a = 0.04, and / = 0.18/ c . In both cases, the full curve corresponds to the 
primary vortex; the dotted curve, to the secondary vortex, (a) Voltage traces V p (t) 
and V s (t). (b) Schematic of time-dependent vortex position (x is the x-coordinate of 
vortex position; a is the lattice constant. 
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Fig. 1 
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Fig. 3 
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